Local electronic properties in nanoscale systems 
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The local electronic structure on nanoscale chains is investigated theoretically. We propose a 
mechanism to explain the even-odd oscillation of length distribution of atom chains. We study 
the spatial peak structure as obtained by scanning tunneling microscopy (STM) constant-current 
topography as a function of the electron-electron interaction, band filling and temperature. The 
site-dependent magnetic moment is also examined. 
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I. INTRODUCTION 

The recent development of techniques capable of 
constructing structures on the nanoscale, with dimen- 
sions that are intermediate in size between isolated 
atoms/molecules and bulk materials, has opened up nu- 
merous possibilities for constructing new devices. This 
perspective raises the importance about novel and im- 
proved physical, chemical, mechanical, and biological 
properties of nanoscale atomic clusters. In particular, 
a relatively recent focus has been to study the effects of 
the electronic quantum confinement within finite length 
chainsii^i^iii^i^. With the help of scanning tunneling 
microscopy (STM) the formation of quantized chain- 
localized states in the pseudogap of the substrate bulk 
band in chains at semiconductor or insulator surfaces has 
been observed^. 

The investigation of the electronic properties of 
nanoscale atomic clusters is particularly interesting. The 
electronic properties suffer strong effects due to the finite 
length of chains. The breaking of translational symme- 
try creates electronic end states. The electronic quantum 
confinement in finite length segments changes the density 
of states within the chains. The conductivity along the 
chain shows spatial variations in the electronic density of 
statesi^^iS. The direct observation of the local electronic 
structure on a nanoscale atomic chain by STM stresses 
the importance of exploring these effects theoretically. 

We can observe that while some STM measure- 
ments can be explained on the basis of non-interacting 
electronsii^i^, others instead suggest strong interaction 
between electrons^ indicating the possibility of a descrip- 
tion using the Hubbard model. The significance of the 
results for extremely small systems without a systematic 
analysis of the observed effects for larger systems would 
not have been appropriate some years ago, when the in- 
vestigation of new nanoscale devices was known only as 
a theoretical possibility. Since the experimental skills 
to manipulate the nanostructures, this study has become 
important, not only to obtain new insight regarding prop- 
erties of theoretical models, but certainly also to better 
understand physical properties in this class of systems. 



In this paper we study theoretically the local behavior 
of the physical quantities of open chains on the nanoscale. 
We report on site-dependent properties of the electronic 
quantum confinement within finite length chains. This 
theory allows for an interpretation of experimental results 
on finite-size effects. We studied two different aspects 
of these systems. First, considering only the electronic 
structure, we introduce a model that captures qualita- 
tively the experimental chain lengths distribution. For 
this case, we analyzed the internal structure of atom 
chains by examining the site-dependent occupation num- 
ber along the axis of the chain for one electron. Sec- 
ondly, we studied the electron-electron interaction by us- 
ing the Hubbard model to understand the main effects of 
the electronic quantum confinement within finite length 
chains. 

The organization of this paper is as follows. Our inves- 
tigation of the chain lengths distribution using the tight 
binding approach is presented in Sec. II. The results 
for local occupation numbers and for the site-dependent 
magnetic moment as a function of the electron-electron 
interaction, band filling and temperature on electronic 
quantum confinement within finite length chains are pre- 
sented in Sec. HI. Our conclusions are presented in Sec. 
IV. 



II. TIGHT BINDING APPROACH 
A. Introduction 

Fitting of experimental wave vectors of nanoscale 
atomic clusters have shown a one-dimensional (ID) free 
electron band dispersion relation^'^. The ID quantum 
potential well has been proposed to describe the confine- 
ment of the electrons and the local conductivity deter- 
mined by the superposition of wave functionsii^i^. In 
contrast to a superposition of wave functions of the ID 
potential wellii^i^, the experimental results are described 
here by the site-dependent occupation number along of 
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FIG. 2: Experimental (full circle), extracted from Ref.— and 
FIG. 1: Experimental (full line), extracted from Ref.^ and theoretical (open circle) distribution of chain lengths for 
theoretical (dotted line) eigenenergies for chains of (a) three, t/{kBT) = 0.5 and (a) low defect density (p — 0.0625) and 
(b) five, (c) seven, (d) nine and (e) fifteen atoms, (e) The (b) medium defect density (p = 0.125). 
electron hopping strength t as function of the chain length. 



the chain. We have used the tight binding Hamiltonian 
T-Cq -t^{cl^Ci+ia + he), (1) 

ia 

where cl^{cia) are the creation (annihilation) operators 
for the electrons of spin a at site i and t is the near- 
est neighbor hopping integral representing the overlap of 
electron wave functions. To gain a quantitative picture 
of the on site-dependence on a chain of N atoms, we have 
computed the local occupation number 

= (*„.|n,|^'„)/7V, (2) 

where n, = hi^ + riij, fiia = c\^Cia and Y^m) is the mth 
eigenvector with energy of Hq. 



As illustration of our approach, we analyze the experi- 
mental data of the eigenenergies for chains of three, five, 
seven, nine and fifteen atoms extracted of Refj^. We have 
fitted the data taking the experimental ground-state en- 
ergies and considering in the tight binding Hamiltonian 
of Eq. 1, conveniently to obtain the experimental band- 
widths. The theoretical and experimental energy spectra 
are shown in Figs. [T^-e. The experimental data are re- 
produced. The electron hopping strength as function of 
the chain length can be observed in Fig. [Tf. Note that, 
as the chain length increases the electronic hopping de- 
creases. The change of hopping amplitudes should be 
caused by changes of the underlying atom chains. 

In the following will investigate the question of the 
length distribution of atom chains. 
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B. Length distribution 

Recent results on the distribution of chain segments of 
gold deposited on a Si surface reveal a relation between 
the chain length and the cohesive energjs^. The size dis- 
tribution is characterized by a strong peak for a length 
of one atom and even-odd oscillations, where even chain 
lengths are favored over odd lengths. We have obtained a 
direct relation between this experimental result and the 
local electronic density of the tight binding model on a 
finite open chain. We find that on some sites the prob- 
ability of finding electrons in the excited energy levels 
is zero. We also observe that these sites correspond to 
positions that divide the chain in m subchains of equal 
length. The delocalized electronic structure determines 
the stability of the atom chain. The cohesive energies 
are a consequence of the local electronic structur o^'^" . 
In particular, we observe that the sites with zero prob- 
ability to find electrons on odd chains are not binding, 
and result in a larger probability for breaking chains. To 
address this problem we propose a simple model. First, 
we consider the distribution of chain lengths for a com- 
pletely random distribution of defects to be p{l — p)^ , 
where p is the defect density^. Since chains are fabricated 
at temperature T ~ IOOO°C and considering i ~ 1 eV, 
it is relevant to evaluate the canonical average for any 
quantity x as 



(x) = Xiexp{~EjkBT)/Z, 



(3) 



where Z is the partition function. Next, we assume that 
there are quantum states breaking a chain of length m 
into g pieces of length N. It follows that g = gm,N — 
{m + 1)/{N + 1) and we can write the distribution of 
chain lengths as 



P(iV) = A^p(l-pr((?(m,iV)), 



(4) 



m— 1 



where ^ is a normalization constant, qi{m, N) = 
gk.N^kN^ k = k{m,i) is a non-trivial function obtained 
numerically, and 5kN is the Kronecker's delta. 

Fig. [2] shows the experimental data, extracted from 
Ref.^, and our theoretical chain lengths distribution for 
low and medium defect densities. The even-odd oscilla- 
tions are clearly visible. Notably, the present model cap- 
tures qualitatively the experimental chain lengths distri- 
bution, considering only the electronic structure effects. 
The quantitative disagreement between theory and ex- 
periment is due to the coupling of the electronic structure 
with the substrate. This coupling was observed in Ref- 
using photocmission measurements of the electronic scat- 
tering vectors at the Fermi surface of the surface states. 

Furthermore, a mechanism for the change in chain 
length obtained by comparing STM images of the same 
sample region taken at different voltage^ is also confirmed 
by our results. Taking different voltages in STM, the 
system can fall into an energy level that singles out to a 
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FIG. 3: Values of no(i) for a chain of 8 sites in the quarter- 
filled band (a) versus U /t for temperature T = and (b) 
versus T for U/t = 8. 



position, which divides the chain in subchains. An ad- 
ditional nearest-neighbor interaction, reflecting the end 
states in the chains, was used to describe this effect in 
Ref.^ and has been already observed in other STM to- 
pography dataiiiS.. 



III. ELECTRON-ELECTRON INTERACTION 

In addition to the tight binding model, containing one 
electron, it is possible to vary the band filling of chain 
structures^. A strong electronic interaction was found 
for small chains of Mn on CuN, and the STM data were 
consistent with the Heisenberg model results^. In this 
case, the effect of electron-electron interaction is partic- 
ularly important. 

The possibility of explaining some STM measurements 
by the tight binding model^'^'^ and others measurements 
by the Heisenberg model^ indicates that we can use a 
description using the Hubbard model. The ID Hubbard 
model is defined by the Hamiltonian 



(5) 



where U is the on-site Coulomb (electron-electron) inter- 
action. 

Solutions for the ID Hubbard model are known since 
the sixtiesi^. The local behavior of the physical quanti- 
ties is however not completely known. Only partial in- 
formations are availabl o^'^'^^i^^ . For example, the spatial 
dependence of the occupation number and of the magne- 
tization was calculated for the ID Hubbard model with 
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open boundary conditions at zero temperature combin- 
ing numerical computations from density-matrix renor- 
malization group and Bethe-Ansatz methoda^^. It is im- 
portant to observe that Friedel oscillations for the density 
and the magnetization in open Hubbard chains have been 
obtained in previous works (see^^ and references therein), 
however, the focus of these studies has been to analyze 
the physical properties of systems having macroscopic 
size. 

We perform exact calculations of physical quantities 
of ID systems having N atoms described by the Hub- 
bard model. We use the standard direct diagonalization 
method^^'^*' and impose open boundary conditions in or- 
der to break the translational symmetry. 

We have found in a half-filled band for U — that 
the occupation number of the ground state is site- 
independent For U > this result is valid for all states. 
This is relevant in order to evaluate the thermodynamic 
quantity (n(i)). We obtain in a half-filled band that 
{n{i)) = 1/N for all sites i, temperatures T and cou- 
plings U > 0. This is a valuable result because the sites 
are not equivalent due to the use of open boundary con- 
ditions. For other than half-filled bands, we observe a 
site-dependence of n,„(i). For a quarter-filled Hubbard 
band we have found a rich dependence on U , m, i and T. 
Fig. [3^ presents the topography of no{i) versus U/t for 
a chain of 8 sites in the quarter-filled band for the tem- 
perature T = 0. At low values of U/t, the intermediate 
atoms have a greater value than average and the end and 
central atoms have smaller values. Here, the probability 
of finding electrons on the intermediate atoms is higher. 
The site-dependence decreases if the electron-electron in- 
teraction increases. For U/t — > oo the site-dependence 
disappears and for all i we find no{i) = 1/8. We can 
alternatively fix U/t and vary the temperature in order 
to cover the other energy states. Fig. [3]d shows this for 
U/t = 8 (other values of U/t give similar results). The 
effect of increasing the temperature is equivalent to an 
increase of U/t. Large temperature destroys the site- 
dependence of the charges at the chain. It is important 
to comment that while this statement appears true for 
the local occupation number, there are other quantities, 
such as intersite spin correlations, that may not follow 
this rule. 

For the magnetic properties of finite chains, first- 
principle calculations have shown that for small chains 
the spin moment depends on position and cluster length. 
For example, for Co chains on Pt(lll), the spin mo- 
ments at the end atoms are higher than those at the 
central atoms and the spin moment of the central atom 
was found to decrease if the chain size increases^ Results 
of the present study are illustrated in Fig. 31 We have 
studied the local spin number 

Sm{^) = 1^™), (6) 

where = (n^i — nii)/2 at site i. 

The local spin number So{i) versus the position i for 
a chain of 7 sites, for the half-filled band and temper- 



ature T = 0, can be observed for different total spins 
S = J2tSoii) and U/t = 8. For S = 1/2 the alternate 
sign of So{i) for odd and even position gives evidence for 
an antiferromagnetic structure. The cases S > 1/2 have 
a ferromagnetic structure, since all S'o(i) have the same 
sign. Site-independence is found for S = 7/2. On the 
other hand, when S is between the minimal and max- 
imal values we find a site-dependence for So{i), whose 
amplitude still alternates but not sufficiently to change 
its sign. For S = 5/2 we obtain that So{i) is higher at 
the end atoms than at the central atoms as has been 
found for Co chains on Pt(lll)^. An important as- 
pect is the influence of the band filling and the clus- 
ter length. In particular, we have found that for the 
case of even total electron number = X^i '^m(*)) the 
Sm{i) is site-independent in all m states. {S{i)) = for 
all sites i, temperatures T and couplings U. Consider- 
ing N even, using that {n{i)) = 1/N in the half-filled 
band and {S{i)) = 0, we have that the canonical ensem- 
ble averages are site-independent, like in mean- field, and 
(n|(«)) — {ni{i)) = 1/{2N) indicating that the open fi- 
nite even chain is paramagnetic in the half-filled band. It 
is important to mention that we find site-independence 
although the system has no translational symmetry. The 
situation for odd finite chains is different. In this case 
(n(i)) = 1/A^, but {S{i)) has a complex site-dependence. 
For the ground state, aX U — Q we obtain a structure 
in which 5o(j) = 1/{N + 1) and on the nearest-neighbor 
sites S'o(i) = {i2^So{i) = 5g = 1/2). At t/ > we 
observe an antiferromagnetic structure. 

Finally, we explore the dependence of the site- 
dependent spin on the coupling U . Fig. shows 5*0 («) 
for the quarter-filled band versus i for a chain of 6 atoms 
and typical U/t values. The ground state has S = 1/2. 
Here, we show the case 5^ = 1/2. For low U/t all SQ{i) 
have the same sign favoring a ferromagnetic order and 
the higher Sq{i) values are at intermediate atoms and 
the smaller at the central atoms. Increasing U/t the cen- 
tral atoms assume negative values for the spin and the 
higher S'o(«) shift to the end atoms following a ttiiTT 
magnetic stripe structure. 

The control and manipulation of the spin rather than 
the charge of electrons^ in STM could show the very 
rich behavior of the local magnetic properties presented 
here. It would be of interest to explore the possibility of 
obtaining experimental results in the context of magnetic 
force microscopy (MFM) devicea^^. This would generate 
images that are sensitive to the spin polarization of the 
electrons. The site-dependent magnetic moment (spin) 
shown here could be confirmed by these experiments. 



IV. CONCLUSION 

We have shown that the spatial peak structure of STM 
constant-current topography can be explained in terms of 
the site-dependent occupation number. We interpret the 
even-odd oscillation of the length distribution of atom 
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chains. We also have studied the spatial structure of the 
electronic charge and spin on finite chains, involving band 
filling, temperature and electron-electron interaction and 
we examined the site-dependent magnetic moment (spin) 
indicating properties that may be confirmed by the MFM 
device. This work provides a simple novel theoretical ex- 
planation for the electronic confinement in atomic chains 
and opens up a way for a more detailed analysis of the 
recent STM measurements. 
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